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We investigate effects of ordinary nonmagnetic disorder in the bulk of a superconductor on magnetic adatom- 
induced Shiba states and on the proximity-induced superconductivity in a nanowire that is tunnel coupled to the 
hulk superconductor. Within the formalism of self-consistent Bom approximation, we show that, contrary to 
widespread belief, the proximity-induced topological superconductivity can be adversely affected by the hulk 
superconducting disorder even in the absence of any disorder in the nanowire (or the superconductor-nanowire 
interface) when the proximity tunnel coupling is strong. In particular, bulk disorder can effectively randomize 
the Shiba-state energies. In the case of a proximate semiconductor nanowire, we numerically compute the 
dependence of the effective disorder and pairing gap induced on the wire as a function of the semiconductor- 
superconductor tunnel coupling. We find that the scaling exponent of the induced disorder with respect to 
coupling is always larger than that of the induced gap, implying that at weak coupling, the proximity-induced 
pairing gap dominates, whereas at strong coupling, the induced disorder dominates. These findings bring out the 
importance of improving the quality of the bulk superconductor itself (in addition to the quality of the nanowire 
and the interface) in the experimental search for solid-state Majorana fermions in proximity-coupled hybrid 
structures and, in particular, points out the pitfall of pursuing strong coupling between the semiconductor and 
the superconductor in a goal toward having a large proximity gap. In particular, our work establishes that the 
bulk superconductor in strongly coupled hybrid systems for Majorana studies must be in the ultraclean limit, 
since otherwise the bulk disorder is likely to completely suppress all induced topological superconductivity 
effects. 
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I. INTRODUCTION 

Majorana fermions in solid-state systems [1-5] obey non- 
Abelian braiding statistics [6] and are a promising platform 
for topological quantum computation [7]. A feasible route 
towards realizing them utilizes a hybrid structure involving 
the proximitization of a semiconductor (SM) with a bulk s- 
wave superconductor (SC) [8-14]. With the appropriate com¬ 
bination of spin-orbit coupling (SOC), Zeeman spin splitting, 
and SC pairing terms, the proximitized system becomes topo¬ 
logical (i.e., an effectively spinless p-wave superconductor) 
and localized Majorana fermions emerge at the ends of one¬ 
dimensional (nanowire) systems or at the vortices of two- 
dimensional systems. Their existence can then be probed 
by conductance measurements as quantized zero-bias peaks 
of height 2e^/h at zero temperature associated with the per¬ 
fect Andreev reflection induced by the Majorana zero energy 
modes [10, 15-19]. Shortly after the theoretical proposals [8- 
13, 20] were put forward, several experimental groups imple¬ 
mented different variants of the proposed Majorana experi¬ 
ment using nanowires in proximity to bulk superconductors 
[21-27]. Although the initial data reporting zero-bias tunnel¬ 
ing conductance peaks in nanowires (albeit with conductance 
values below the theoretically predicted 2e^/h quantized con¬ 
ductance) are encouraging, more theoretical and experimental 
work still needs to be done in order to distinguish signatures of 
Majorana fermions from those from other possible nontopo- 
logical mechanisms as have been discussed in the literature 
[28-35]. 

The current work is on the deleterious effect of disorder 
on the proximity-induced topological superconductivity in the 
hybrid system of experimental interest. The topological super¬ 


conductivity induced in the nanowire, arising from a combi¬ 
nation of s-wave superconductivity, spin splitting, and spin- 
orbit coupling, is essentially equivalent to a type of an ef¬ 
fectively spinless p-wave superconductivity [36] with triplet 
spin correlations [37] which is not immune to ordinary non¬ 
magnetic disorder in the environment, unlike regular i-wave 
spin-singlet superconductors which are protected against non¬ 
magnetic disorder by virtue of the Anderson theorem. There 
have therefore been many theoretical and numerical studies 
[28-30, 38-54] of the effects of disorder on the topological su¬ 
perconductivity in this context, going back to almost 15 years 
ago [28]. It may appear that another theoretical study of dis¬ 
order effects in this context would be redundant, but as we 
explain below, this is not the case here. The specific question 
regarding disorder effects (in the bulk superconductor itself) 
addressed in this paper has only been discussed a few times 
in the literature, with the first paper [55] coming to an erro¬ 
neous conclusion which was subsequently corrected [56, 57]. 
The conclusion we reach in our current work is of great im¬ 
portance in choosing the proper materials for the hybrid struc¬ 
tures manifesting topological superconductivity and Majorana 
fermions. 

In considering the effects of disorder in these hybrid SM- 
SC [21] (or FM-SC where FM stands for the ferromagnetic 
adatoms as in Ref. [27]) systems, one should distinguish be¬ 
tween disorder in the SM and that in the bulk SC. Disorder in 
the SM has been extensively studied [30, 39, 40, 44, 46, 47], 
with the main conclusion being that the topological gap is de¬ 
stroyed by this type of disorder when the mean free path is 
comparable with or smaller than the induced coherence length 
in the semiconductor. It has been much emphasized in the lit¬ 
erature [44] that the topological nanowire must be in the bal- 
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listic limit with the carrier mean free path being much larger 
than the proximity-induced coherence length in the nanowire 
for the manifestation of the Majorana zero modes, a condi¬ 
tion which is likely (unlikely) to be satisfied in the semicon¬ 
ductor [21] (ferromagnetic [27]) nanowires. It has also been 
emphasized [58] that the applicable disorder at the nanowire- 
superconductor interface must be low for the induced proxim¬ 
ity superconductivity to manifest a hard gap as has recently 
been reportedly accomplished in the InAs-Al epitaxial core 
shell nanowire hybrid structures [59]. On the other hand, 
disorder in the bulk SC has received relatively little atten¬ 
tion [14, 55-57], with the focus mainly on the limit where 
the coupling between the two materials is small (i.e. the 
weak-coupling limit where the SM-SC tunneling amplitude 
is small). In this limit, it has been found [55, 57] that the dis¬ 
order in the bulk SC hardly affects the superconducting gap in 
the topological system, making it possible to use disordered 
or dirty SCs in experiments [21-25]. One consensus in the 
community regarding disorder effects seems to be that disor¬ 
der in the nanowire itself (superconductor itself) is important 
(unimportant) with respect to the manifestation of proximity- 
induced topological superconductivity and Majorana fermions 
in the hybrid system. The current work directly challenges 
this consensus, showing that the disorder in the bulk supercon¬ 
ductor may very well be important for the proximity-induced 
topological superconductivity, particularly in the limit where 
the superconductor and the nanowire are strongly tunnel cou¬ 
pled. In particular, the bulk superconductor should be in the 
clean limit, with its elastic mean free path being much larger 
than the superconducting coherence length for optimal in¬ 
duced topological superconducting order in the hybrid struc¬ 
tures. Our current work indicates that having a clean su¬ 
perconductor with a very long mean free path is an absolute 
necessary condition for the realization of a large proximity- 
induced topological superconducting gap hosting Majorana 
fermions in the SM-SC and FM-SC hybrid systems. 

Two recent developments prompted us to revisit the issue 
of bulk disorder in the superconductor. First, a new class 
of proposals [27, 60-73], utilizing Shiba states induced by 
magnetic adatoms on SCs to generate Majorana fermions, has 
emerged. This platform for using ferromagnetic adatoms on 
a superconducting substrate as the topological FM-SC hybrid 
system is, in some sense, the large spin splitting limit [34] of 
the SM-SC hybrid structure with the spin-splitting arising in¬ 
trinsically from exchange effects in the ferromagnet, rather 
than from a Zeeman splitting induced by an external mag¬ 
netic field as in the SM-SC hybrid system. This ferromagnet- 
superconductor hybrid system can therefore be effectively de¬ 
scribed by a Hamiltonian which is the same as that of the SM- 
SC heterostructure, with a crucial difference that the tunnel 
coupling between the adatoms and the SC (as well as the spin 
splitting in the adatom chain) is much larger than the corre¬ 
sponding term in the SM-SC system [35], rendering the pre¬ 
vious perturbative treatment of disorder in SC inapplicable. 
Second, in Ref. [59], the SM-SC structure has been grown 
epitaxially, which drastically improved the quality of the inter¬ 
face between the two materials. A hard proximity-induced su¬ 
perconducting gap is then observed on the SM, resolving the 


“soft-gap” issue that previous experiments found [58]. The 
size of the gap on SM is comparable to that on the SC, in¬ 
dicating strong coupling between the two materials [74]. In 
this limit, however, it is unclear whether disorder in the bulk 
SC can significantly degrade the gap on the SM, especially 
when a magnetic field is applied on the SM to create a Zee- 
man spin splitting necessary for producing the topological su¬ 
perconductivity in the SM wire. 

In both of these experiments [27, 59], the strong tunnel cou¬ 
pling between the SM (or magnetic adatoms) and the SC ne¬ 
cessitates the re-examination of the issue of disorder in the SC, 
as previous treatments of this problem were valid only when 
the coupling is small [55-57], which is not the case in these 
two new systems. In this paper, we investigate the effects of 
disorder in a SC on the spectral properties of a proximate SM 
in the SM-SC system and on the Shiba states in the FM-SC 
system. The disorder problems for the two hybrid structures 
(i.e., disorder effects on the Shiba states in the ferromagnetic 
adatom chain and on the SC/SM nanowire) are somewhat dif¬ 
ferent, and we therefore study the two systems (FM/SC and 
SM/SC) separately so that our work applies to both experi¬ 
mental systems, although our main emphasis in the current 
work is on the semiconductor-based Majorana hybrid systems 
since the experimental situation is better understood in such 
semiconductor-nanowire structures. The formalism we adopt 
is the self-consistent Born approximation, which is valid in the 
limit of weak impurity scattering (specifically kfl » 1, where 
kf and I are, respectively, the Fermi wave number and the 
disorder-induced transport mean free path in the bulk SC- this 
condition is well satisfied in the bulk superconductors used 
in the Majorana hybrid structures with the clean/dirty bulk 
superconductors being defined by whether 1 ^ or 1 

respectively, where f is the SC coherence length). We extract 
the density of states (DOS) of the topological systems via their 
dynamic Green functions which contain the self-energy due to 
ensemble-averaged disorder in the bulk of the SC. Throughout 
the paper, we shall assume the SM itself as well as the SM- 
SC interface is clean and only consider disorder in the bulk 
superconductor. 

Investigation of whether a strong tunnel coupling to the su¬ 
perconductor in hybrid systems, e.g. the inherent strong cou¬ 
pling in metal-on-metal ferromagnetic adatom systems or the 
strong coupling in epitaxial semiconductor-superconductor 
systems, could lead to the disorder in the bulk superconduc¬ 
tor (quite apart from the disorder in the adatoms or the semi¬ 
conductor itself) becoming relevant for the topological super¬ 
conducting (and consequently Majorana fermion) properties 
is the goal of the current theoretical work. We ignore all dis¬ 
order in the SM nanowire (or the FM chain) itself since it has 
already been studied extensively elsewhere and is well under¬ 
stood as being detrimental to topological properties. 

To qualitatively understand how the bulk superconducting 
disorder might be relevant for the proximity-induced super¬ 
conductivity in the strong tunnel coupling (which we often re¬ 
fer to simply as “strong coupling” in this paper) limit, we refer 
to the simple (and approximate) formula for the proximity- 
induced superconducting gap in the hybrid system derived in 
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Refs. [10, 44, 74] and used extensively: 




TA 

r +A 


( 1 ) 


where F is the effective coupling, A is the bulk gap in the 
superconductor, and A,,, is the induced proximity gap in the 
nanowire. It has been much emphasized that in order to ob¬ 
tain a large induced gap, one must have F » A so that A„ ~ A, 
which is obviously the maximum possible value of the in¬ 
duced gap (as achieved presumably in the epitaxial InAs-Al 
core-shell nanowire systems [59]). In the opposite limit of 
very weak coupling, F A, one gets A„, ~ F with a very 
small induced gap « A (with consequently an even smaller 
topological gap since the topological gap is bounded from 
above by A^,). Let us now imagine an extremely large tun¬ 
nel coupling (e.g., F going to infinity) where there is then 
no discernible difference between the superconductor and the 
nanowire so that A„ = A applies, and hence the nanowire has 
essentially become a part of the bulk superconductor as far as 
superconducting properties go. In such a situation, the bulk 
disorder in the superconductor is now a part of the disorder 
in the nanowire since from the perspective of superconduc¬ 
tivity, these two have become one monolithic system. Now, 
if we turn on spin-orbit coupling and spin splitting so as to 
convert A„ into a topological superconducting gap, then the 
disorder existing in the bulk superconductor must necessarily 
suppress the effectively triplet spinless topological supercon¬ 
ductivity since it is not protected by any Anderson theorem 
(as time reversal invariance is explicitly broken). We note that 
this argument does not apply in the weak tunneling (F « A) 
limit where the two parts of the hybrid system (the bulk su¬ 
perconductor and the nanowire) are distinct, and indeed it has 
been explicitly shown [57] that in the weak tunneling limit 
(F going to zero), the bulk superconducting disorder does not 
suppress the topological superconductivity, but of course the 
topological gap is very small in this limit (< F) anyway. Al¬ 
though this physically motivated qualitative argument is not 
a proof by any means, the argument demonstrates that the 
strong- and weak-coupling situations could be fundamentally 
different with respect to disorder effects coming from the bulk 
superconductor, and a careful investigation is necessary to see 
whether or not the strong-coupling situation is benign with 
respect to the bulk disorder effects. While this problem is of 
considerable intrinsic interest itself in the context of the theory 
of superconducting proximity effect, the current experimental 
push (e.g., InAs-Al epitaxial hybrid system) to produce a hard 
induced gap makes our work timely in the study of Majorana 
fermions in solid-state systems. 

The paper is organized as follows. In Sec. II, we consider 
the ferromagnet-SC hybrid system in the Shiba limit, where 
the coupling between the magnetic adatoms and the SC is 
much stronger than the interatomic coupling. We find that 
in this strong-coupling regime, disorder in the bulk SC has 
strong effects on the location of the Shiba energy in the bulk 
SC gap. In Sec. Ill, we consider the effects of bulk supercon¬ 
ducting disorder on the SM-SC heterostructure, in both the 
weak- and strong-coupling limits. We show that our results in 
the weak-coupling limit agree with previous works [55, 57], 


and highlight features specific to the strong-coupling limit, 
where, in contrast to the weak-coupling limit, nonmagnetic 
elastic disorder in the bulk superconductor invariably strongly 
degrades the proximity-induced topological SM superconduc¬ 
tivity. We conclude in Sec. IV with a summary and with a 
brief discussion of the far-reaching implications of our find¬ 
ings for the future design of hybrid structures hosting Majo¬ 
rana fermions. 


II. FERROMAGNETIC ADATOM-INDUCED SHIBA 
STATES IN A DISORDERED SUPERCONDUCTOR 

We first consider a disordered s-wave SC strongly coupled 
with a magnetic impurity, described by 

kcr k 

-J ^ a-al (r = 0) Oc- (r = 0) (2) 

(T 

+ [+^ drUdis (r) ^ a'l (r) (r), 

cr 

where Uka- annihilates an electron with momentum k and spin 
cr. In the first line, is the normal-state dispersion and A is 
the i-wave pairing term of the SC. In the second line, J char¬ 
acterizes the strength of the magnetic impurity, located at the 
origin (r = 0), which induces a local Zeeman term in the SC. 
The prefactor cr = +1 corresponds to spin-up/down respec¬ 
tively. In the third line, t/dis (r) represents nonmagnetic elas¬ 
tic disorder present in the SC (which leads to a finite transport 
mean free path / in the bulk SC in its normal state). Below, we 
investigate effects of the magnetic term (7) on the DOS of the 
system, for both clean (t/dis = 0) and dirty (t/dis 7= 0) SCs by 
generalizing the original Yu-Shiba-Rusinov theory [75-77] to 
include static nonmagnetic elastic disorder t/dis in the SC. 


A. Clean Superconductor 


We first briefly review the theory of Shiba states [75-77] in 
the absence of disorder (t/dis = 0) to set a context and to fix 
the terminology. In frequency (cu)-momentum (k) space, the 
Green function for the system is 

^kk' ~ ^k ^ fcfc'G^y, fj; 


where (ai) = Green function for a clean 

sc with Tfj the Pauli matrices acting on the Nambu-Gorkov 

space of . The superscript (1) in Eq. (3) indicates 

the presence of one magnetic impurity, but without disorder in 
the system. The effect of the magnetic term J is captured by 
the T matrix in the second term of Eq. (3), which is given by 


Tkk' 





J 

V 


^1 - Jttvo 


OJTo + ATj 

- A^ 


J 

V’ 


(4) 
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Figure 1. (Color online) The LDOS at the position of the magnetic 
impurity in a clean (black solid lines) or disordered (colored lines) 
SC. The Shiba-state energies are tuned to (a) eo = 0 and (b) eo = 
0.4A. A broadening of magnitude O.OOIA is used to smear out the 6- 
functions for depicting the results. The elastic disorder in the system 
is quantified by the mean free path I, which is given in the units of 
the clean limit coherence length ^ of the SC, as shown on the right. 


where V is the volume of the system and vq is the normal-state 
DOS at the Fermi level. The pole of T in the subgap regime 
indicates the presence of a bound state, called the Shiba state 
[75-77], with the energy given by 


eo = sgn7 


1 - (Jnvof 
1 + (jTTVof 


A. 


(5) 


The local density of states (LDOS) at the position of the Shiba 
state is given by v(r = 0) = Zifcfc' Pig- 1’ the 

black lines show the the LDOS at a Shiba state with ener¬ 
gies eo = 0 and eo = 0.4A. The divergence of the LDOS at 
w = eo indicates that the Shiba states have well-defined ener¬ 
gies. We mention the obvious fact that the Shiba-state energy 
eo is tuned by appropriately tuning the magnetic coupling J, 
and in a given experimental setup, J would typically be fixed 
producing a Shiba energy according to Eq. (5) above. Re¬ 
sults corresponding to two situations with eo = 0 and 0.4A are 
shown in Figs. 1 and 2. 


B. Disordered Superconductor 

We now investigate the effects of ensemble-averaged dis¬ 
order in the bulk SC (t/dis) on the Shiba-state energy. To this 
end, we assume that the random quenched nonmagnetic impu¬ 
rities in the SC have a concentration of Wimp, and each impurity 
has a scattering potential of the form t/, (r) - U6 (r - ri), 
where r, is the position of the impurity. In the self- 
consistent Born approximation, the Green function is written 


as 

Gkk' = [(C^fc,) - Sfcfc/ 


( 6 ) 


where is given by Eq. (3) and the inversion here is oper¬ 
ated on the k - k' matrix space. The disorder-induced self¬ 
energy Sfcfc/ is given (in the self-consistent Born approxima¬ 
tion) by 


"^kk' — ^impG ^ ^ 

P 

T“' 1 V 


(7) 

( 8 ) 


where t = (27rninipG^vo) is the disorder scattering time. 
Note that due to the lack of translational invariance, the self¬ 
energy due to disorder, in its exact form of Eq. (7), is nondi¬ 
agonal in fe-fe'. It is easy to check that Eq. (7) reduces to the 
conventional form for translationally invariant systems [78] if 
Gkk' is proportional to 6kk'- 

In reaching Eq. (8), we observed that that has the 
highest weight when |fc| = |fe'| = \ where is the Fermi 

momentum of the system in its normal state. Therefore, in the 
summation over p in Eq. (7), the summand has appreciable 
weights only when \p + k\ - |p + fc'| = In the general 
case where k + k', this condition is satisfied only for a one¬ 
dimensional manifold of p, but when k = k', the condition 
reduces to \p + k\ - and is satisfied by a two-dimensional 
manifold of p. Thus we see that I,kk' has most of its weights 
at A; = k', allowing us to approximate it by Eq. (8). This can 
be understood physically as the ensemble-averaged disorder 
should not introduce further translational-symmetry breaking 
and hence I-kk' is diagonal in k. 

For each value of energy (w), Eqs. (6) and (8) are iter¬ 
ated numerically until convergence. In the evaluation of the 
momentum integrals, we use the approximation ^ ^ 

Vq d^k- The ^k integral is discretized on a grid with 
10^ points distributed in a way such that sampled uni¬ 

formly over the interval |^10 ij. The iteration converges af¬ 
ter a few cycles for most values of a> except for those in the 
vicinity of the Shiba-state energy (eo), which requires a few 
hundred of iteration cycles. 

The blue dashed line in Fig. 1 shows the LDOS at the mag¬ 
netic impurity when the SC is disordered with a mean free 
path of / = VfT - lOf, where Vf and ^ = Vf/A are the Fermi 
velocity and the coherence length of the SC, respectively. The 
results for several other values of / are also presented. We ob¬ 
serve that the 6 peak at a>/A = eo associated with the Shiba 
state is now broadened to a dome by the disorder in the SC. 
This can be understood as follows: the continuum states of the 
SC are different for each realization of disorder. The scattering 
phase shift due to the magnetic impurity therefore varies from 
one realization to another, leading to an effective disorder- 
induced fluctuation in J and eo [c.f. Eq. (5)]. Ensemble av¬ 
eraging the LDOS over disorder leads to a domelike shape 
centered around the energy of the Shiba state in a clean SC 
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Figure 2. (Color online) (a) The normalized fluctuation SJjJ as a 
function of the disorder strength f// and the Shiba energy eo for clean 
SC. In the white regions, the subgap dome of LDOS joins the contin¬ 
uum modes (|m| > A) which makes 6J ill defined, (b) The width of 
the subgap dome as a function of disorder strength for eo = 0 (solid 
line) and eo = 0.4A (dashed line). 


with the dome in Fig. 1 reflecting the “spreading” in the effec¬ 
tive Shiba energy due to disorder—it is clear that the clean SC 
limit with f ^ is necessary for the system to have a sharp 
Shiba energy. 

We plot in Fig. 2(a) the effective fluctuation in J, which is 
defined as the standard deviation 6J of the distribution in J 
that would result in Shiba-state energies distributed according 
to the subgap LDOS obtained from the self-consistent Born 
approximation (e.g., the dome in Fig. 1). Fixing the Shiba- 
state energy in the clean limit at eo = 0 and 0.4A (as in Fig. 1), 
we plot in Fig. 2(b) the width of the subgap dome of LDOS 
against the strength of disorder as characterized by In 
general, the Shiba subgap state is broadened by disorder, and 
the subgap DOS joins the continuum states when / < with 
the precise critical disorder dependent on the Shiba energy. 
The results presented in Figs. 1 and 2 clearly demonstrate the 
importance of being in the ultraclean-SC limit, i.e., / » for 
obtaining a Shiba state whose energy is close to that in the 
clean limit. 


C. Discussion 

We have investigated the LDOS associated with a single 
magnetic impurity embedded in a disordered SC with ensem¬ 
ble averaging. The result thus obtained is not expected to be 
directly applicable to a real experiment conducted with a sin¬ 


gle impurity since in reality there is only one realization of dis¬ 
order, and therefore this Shiba-state energy should appear as 
a sharp subgap LDOS peak. Our ensemble-averaged results, 
however, reveal that the Shiba-state energy is shifted ran¬ 
domly from sample to sample around its value in the clean-SC 
limit (with the likelihood roughly proportional to the height of 
the subgap dome) because each specific experimental sample 
will have its unique disorder configuration which will differ 
randomly from one sample to another. Qualitatively, when 
the mean free path of the bulk SC is of the same order of 
magnitude as (or shorter than) its SC coherence length, the 
Shiba-state energy could be anywhere within the SC gap. 
This is likely to have implications for a class of recent pro¬ 
posals which utilize the Shiba states induced in a supercon¬ 
ductor by a chain of magnetic adatoms to generate Majorana 
fermions [27, 61, 64, 66, 67, 69, 72, 73]. In general, the spin- 
orbit-coupling strength, the lattice spacing between magnetic 
adatoms, and the Shiba-state energy all need to be fine tuned 
to obtain topological superconductivity in the system [67, 72]. 
However, it is difficult to control the Shiba-state energy even 
if the SC is clean, as it is determined by the strengths of the 
magnetic adatom and its coupling with the SC, both of which 
can hardly be tuned experimentally. Our finding indicates that 
in addition to the difficult task of fine tuning the Shiba energy 
to zero, it is also necessary to use an ultraclean SC with I » f 
so as to ensure that the Shiba-state energy for each individual 
atom remain close to zero. Otherwise, with Shiba-state ener¬ 
gies along the chain being shifted randomly by disorder, the 
system would then be fragmented into segments of topolog¬ 
ical and nontopological regions, an unfavorable situation for 
topological quantum computation. Thus, our current work im¬ 
plies that the fine-tuning problem of creating Majorana modes 
using Shiba states becomes substantially worse in the strong¬ 
coupling situation since the bulk disorder in the supercon¬ 
ductor now randomly shifts the Shiba-state energy, leading to 
strong sample-to-sample variations. We mention here that the 
typical Shiba-induced Majorana system is a metal-on-metal 
system (i.e., a ferromagnetic metal chain on a superconduct¬ 
ing metal) where the tunnel coupling is large (~eV), and the 
typical bulk disorder scale in the superconductor (~meV) is 
much larger than the typical induced gap (~ 0.1 meV), leading 
to an intrinsically unfavorable theoretical situation for the ex¬ 
istence of Majorana modes by virtue of the fine-tuning prob¬ 
lem. 


III. SEMICONDUCTOR 
NANOWIRE-SUPERCONDUCTOR SYSTEM 


In Sec. II, we looked into the system in the “Shiba limit,” in 
which the interatomic hopping among the magnetic adatoms 
is much weaker than their coupling with the bulk SC. We now 
turn to consider the opposite limit where the system hybridizes 
strongly to form a nanowire. (We mention as an aside that 
the two limiting situations, i.e., the Shiba limit of weak inter¬ 
wire hopping [69] and the nanowire limit of strong interwire 
hopping [70], are not separated by a quantum phase transi¬ 
tion and are two extremes, which are applicable to different 
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physical situations, of the same underlying physics [35, 73].) 
Previous work [57] has indicated that disorder in the bulk SC 
cannot degrade the superconducting gap in the nanowire if 
the SM-SC coupling is weak. It, however, remains unclear 
whether a strong SM-SC coupling could alter this conclusion 
qualitatively (as was already discussed in Sec. I of this paper). 
Therefore, we now investigate this question in depth without 
assuming any weak SM-SC coupling. Note that the effect of 
disorder at the SM-SC interface is a separate issue which has 
been theoretically treated previously [58]. 

In the absence of disorder, the system is described by the 
Hamiltonian H - H„ + H^c + Hj, where is the Hamilto¬ 
nian for the SM wire/SC and Hj is the coupling between the 
two materials. Explicitly, they are given by 


k-<T 

Hsc - ^ ^ XI 

k(r k 

kcr 

where - ju and are the dispersions of the wire 

and the SC, respectively. B and a are the Zeeman and SOC 
terms on the wire, and A is the i-wave pairing term on the SC. 
The subscript a represents the spin species opposite to that of 
cr. The physical properties of the system are captured by its 
Green function, which is given by 


Inverting the matrix in Eq. (10), we get the following set of 
coupled equations: 


Gw = (g™-' - E,.,)^ , 

(12a) 

Sw = T X^^^^ (fc)Tz. 

^ kx 

(12b) 

Gs = (Gf-' - Edis)^’, 

(12c) 

k') y X ■^zG^(gi,q 2 )Tj, 

(12d) 

fc') = Gs6kk' + G.S (fc) f^Tj X 


Gn. (kj) TjGj (fc') 6k^^k>, 

(12e) 


where is the Green function of the SC with the effects of 
disorder incorporated, while Gs incorporates both the effect of 
disorder and that of the coupling to the wire. [A in Eq. (12b) 
is the area required for normalization similar to the volume V 
normalization in the earlier equations.] 

The DOS on the wire can be calculated from G„ by 
v„ («) = ^Trlm I?! ^G„ {kf). To evaluate the momentum 
integrals, we use the following approximations 


1 

A 

1 

V 


E 

k± 

E 


—> ^ 2 D [+1 


(13a) 

(13b) 



Here, the full Green function ^ of the whole system is written 
in the block spinor space of where repre¬ 

sents the Bogoliubov-de Gennes spinors for the SC/wire, re¬ 
spectively. The off-diagonal block matrix T represents the 
tunneling between the two systems. Edis is the self-energy 
originating from the nonmagnetic disorder, which is present 
in the SC only in our model. The explicit forms for these 
terms are 


WTO H- f -I- At 


(11a) 


Gf)(«,fc«)_ 

(w,4“'^) = (wTo - Tz - Bo-:. - ak'f'^o-jcT^) ' , (lib) 


(11c) 


^ X Gsiq'fKq'f), (lid) 


where G^^^ are the unperturbed Green functions of the 
wire/SC. The superscripts (w) / (i) on the momentum vari¬ 
ables indicate that they refer to the wire/SC. The effects of 
disorder in the bulk SC are captured by the self-energy Edis, 
whose expression has been approximated in the same way as 
Eq. (8). 


where V 2 d — vo^lkp is the density of states of a two- 
dimensional system and kp is the Eermi momentum. 

We are primarily interested in the spectral properties of 
the SM wire, which becomes topological and hosts Majorana 
fermions with suitable combinations of SOC, Zeeman and in¬ 
duced SC pairing terms. The influence of the SC on the SM 
wire is capmred by the self-energy term E„. To have a better 
understanding, we expand it at small frequencies as 

(oj) ss (Eo + wEq) + EjcTx + ■.., (14) 

where E^ = ^Tr |^E,^, (0) t^I and Eq = jTrReE]^ (0) are scalar 
numbers, and the symbol ..” represents terms proportional 
to other matrices. These terms renormalize the SOC, Zee- 
man splitting, effective mass, and chemical potential of the 
SM wire, and do not concern us here. By substituting Eq. (14) 
into Eq. (12a), we obtain the following form of Green function 
for the wire: 

G„ = [w -- (Eo -I- wEq) - E.tTj:] , (15) 

= Z-^ ^co + ^ - Z-^9{„ - Z-%Txj , (16) 

where Z - 1 - Eq is the frequency renormalization factor, and 

T„' = -Z-'lmEo (17) 

is the broadening induced by Ej^,. Since t“^ has the same ef¬ 
fect as disorder, we define it as the effective disorder on the 
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SM wire which can be thought of as the proximity-induced 
effective disorder arising in the SM due to the presence of 
the SC. The term Z induces a SC pairing onto the wire, 
but we shall not define it as the pairing term directly since 
Eq. (14) is an expansion at zero frequency and does not cap¬ 
ture the frequency dependence of Rather, we numerically 
compute the spectral gap when the wire has no Zeeman split¬ 
ting {B - 0), and identify this gap as the effective magnitude 
of the pairing term (A„) in the SM wire. 

Before presenting the full numerical results, we review the 
conventional treatment of the problem in the weak-coupling 
limit where E A, with E = nv 2 Dt^ being the tunnel coupling 
strength between the wire and the SC. In this limit, the sec¬ 
ond term of Eq. (12e) can be neglected, and Eq. (12) has the 
analytic solution [79, 80] 




WTO -H + ATx 


- A2 




k 

o) — Atv 


Va2“ 


(18) 

(19) 


where d> and A satisfy to — a> + j — ?== ^nd A = A -H 

2— . The self-energy therefore gives a superconduct- 

ing gap of size E (« A) on the wire. Since ImS„, (w = 0) = 0, 
we see that disorder is not induced on the wire in this limit. 
This is the gist of the weak-coupling result obtained ear¬ 
lier by Lutchyn et al. [57] establishing the immunity of the 
proximity-induced topological superconductivity to any dis¬ 
order in the SC itself, and it is only valid for E « A when the 
induced gap is extremely small (~ E). 

In the strong-coupling limit in which the second term of 
Eq. (12e) is not small and cannot therefore be ignored, both 
superconducting gap and disorder are induced on the wire by 
the disordered SC. We now investigate their dependence on 
the strength of the bulk disorder and of the SM-SC coupling. 




Figure 3. (Color online) Induced superconducting gap on the wire 
(A,,,) against SM-SC tunnel coupling (F) in (a) linear scale and (b) 
logarithmic scale. Different lines correspond to different strengths of 


disorder in the bulk of the SC 

wire are chosen to he 6 = a = 0 and p = A. 


The parameters on the 



A. Superconducting Gap Induced on the Wire 

We first investigate the superconducting gap induced on the 
wire which, in the absence of Zeeman term on the wire, is 
equal to the induced pairing potential A^. For simplicity we 
set B - a = 0 and yU = A on the wire. (We have explic¬ 
itly numerically checked that our results presented here are 
generic, and using different parameter values do not change 
the results at all qualitatively.) We note that the natural en¬ 
ergy scale of the problem is A and the natural length scale is 

The dimensionless parameter quantifying disorder is 

defined as d - where I is the transport 

mean free path of the bulk SC. 

In Fig. 3, we plot the calculated induced gap A^ as a func¬ 
tion of E in both linear and logarithmic scales. We observe 
from the log-log plot that in the weak-coupling and weak dis¬ 
order limit, the induced gap scales linearly with coupling, i.e., 

A„ oc r, for I » and E ^ A. (20) 


Figure 4. (Color online) Calculated dimensionless induced gap 
(A„,/F) plotted against the dimensionless tunnel coupling (A„,/F) for 

various values of the disorder parametrized by li = as 

shown. The parameters on the wire are chosen to be 6 = a = 0 
and p = A. 

This is the same as previous calculations done with a clean 
SC [74]. Perturbative calculations [57] have also shown that 
at weak enough coupling T, disorder in the bulk SC does not 
change the linear scaling of A^, with E. To verify this, in Fig. 4 
we plot Avv/E against E, where we observe that 

lim^ = l, (21) 

r^o E 

irrespective of the strength of disorder in the SC, which is 
identical to the induced gap for a clean SC in the weak- 
coupling limit (E « A), as we pointed out in the discussion 
following Eq. (19). This is an indication that the induced pair¬ 
ing dominates over the induced disorder in the weak-coupling 
limit E « A. In the following section, we develop a more 

















8 




d 

—0.01 
- 0.032 
-0.1 

--- 0.32 

-1 

--- 3.2 
—10 


Figure 5. (Color online) Effective disorder strength on the wire 
as a function of SM-SC coupling (F) in (a) linear scale and (b) loga¬ 
rithmic scale. The parameters on the wire are chosen as a = 0, yu = A, 
and B = 5A. Different lines represent different disorder strength de¬ 
fined hy d = yjk^/l. 


quantitative understanding by computing the scaling behavior 
of disorder with respect to F and comparing it with Eq. (20), 
explicitly demonstrating that the induced pairing dominates at 
small r. 


B. Disorder Induced on the Wire 



Figure 6. (Color online) Effective disorder on the wire against 

disorder in bulk SC (C') in logarithmic scale. The parameters on the 
wire are chosen as or = 0, p = A, and B = 5A. The disorder strengths 


are defined 


hyd= yj^/L 


<si r. However when is comparable with F, the de¬ 
pendence changes to linear. Similarly, from the plot of 
against in Fig. 6, we see that it scales linearly when T' is 
small. In summary, we have 


T 


-1 

w 


Ir''’' 


F, 
< F. 


(23) 


This result could be compared with previous results ob¬ 
tained from perturbative treatments which apply in the limit 
of very small induced disorder [55, 57]. We note that in ad¬ 
dition to recovering the quadratic scaling at weak coupling, 
our result also shows a crossover to linear scaling at interme¬ 
diate coupling strength, which cannot be captured by pertur¬ 
bative approaches. This qualitatively different linear scaling 
behavior of the induced disorder has important consequences 
for the Majorana-carrying SC-SM hybrid nanowire systems as 
discussed below. 


Apart from the proximity-induced superconducting gap, 
the wire also inherits disorder from the SC. In the topologi¬ 
cal phase, the superconducting gap protecting the Majorana 
fermions could be destroyed by the induced disorder if its 
strength becomes comparable to the gap [44]. It is therefore 
important to study the dependence of the induced disorder on 
the coupling strength so as to compare its strength with that of 
the gap. The disorder on the wire has been defined in Eq. (17), 
which when written out in full is 

r-' = Trims,, (cu^O) 

TrReS;, (w = 0) - 4 ^ 

We remark here that could, in general, be defined as a 
frequency-dependent quantity, but for our purpose of investi¬ 
gating its dependence on coupling strength we take only its 
value at zero frequency. Also, from Eqs. (12), we see that 
ImS„, {a>) is nonzero only if the band dispersion of the wire 
crosses the energy oj. Therefore, in order to produce a nonzero 
we choose the Zeeman term on the wire to be B = 5A, 
while keeping a = 0 and // = A, noting that the results depend 
only weakly on the parameters. 

In Fig. 5, we plot against F in linear and logarithmic 
scales. It shows that scales quadratically with F when 


C. Suppression of Topological Gap by Disorder 

For practical reasons, the calculations in Secs. IIIB and 
III A have been performed for wires with and without Zee- 
man terms, respectively. It is natural to ask whether there is 
a system in which both induced SC gap and induced disorder 
are at play. The Majorana nanowires, with nonzero SOC and 
Zeeman terms, are such systems. If its normal-state dispersion 
crosses zero energy, the disorder term as defined by Eq. (17) 
is nonzero. On the other hand, with the SOC term present, the 
induced pairing can produce a spectral gap even in the pres¬ 
ence of the Zeeman term. Therefore, a comparison between 
the scaling behaviors of induced pairing and induced disorder 
with respect to the SC-SM tunnel coupling is in order. 

Comparing Eq. (20) and Eq. (23), we see that the scal¬ 
ing exponent (2) of disorder with respect to coupling 
strength (F) is larger than that (1) of the pairing term (A,). 
This implies that A„ dominates at smaller F, while domi¬ 
nates at larger F. Therefore, increasing the coupling strength 
does not always lead to a larger proximity-induced topolog¬ 
ical gap since the induced effective disorder increases faster. 
Rather, the optimal value of F at which the topological gap 
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Figure 7. (Color online) Superconducting gap on the SM wire (£gap) 
as a function of coupling strength (F) for various values of disorder 


strengths in the bulk SC 



The parameters on the wire 


are chosen to be ^ = A and (a) B - 2A, a - 0.3 



= 2A, 


a = 


(c) B = 3A,a = 0.3 (d) S = 3A, or = 



is maximum is dependent on the parameters of the system, 
and one expects some intermediate value of F, which is nei¬ 
ther too small (so that the intrinsic topological gap itself is not 
too small) nor too large (so that the induced disorder is not 
too strong overwhelming the induced gap), to be the optimal 
choice. Purely on dimensional grounds, the optimal value of 
the tunnel coupling is expected to be F ~ A in the relatively 
clean-SC limit, but with increasing SC disorder, we expect the 
optimal value of F to decrease so as to keep induced disorder 
effects small. 

Figure 7 shows the superconducting gap on the wire (£gap) 
a s a funct ion F for several values of B and a. For 0 < F < 
- the system is in a topological regime since the s- 
wave pairing term is smaller than the Zeeman term [11, 13]. 
When a is nonzero and in the absence of disorder (black line 
in Fig. 7), a topological gap exists and a zero-energy Majo- 
rana fermion is present at each end of the nanowire. When 
F increases to values near the Zeeman energy |B|, the topo¬ 
logical gap shrinks and closes completely at F = 
indicating the onset of the topological phase transition. The 
gap then reopens at F > ^jB^ - where the system is non- 
topological. 

The gap on the wire is suppressed by disorder in the bulk 
SC, but the degree of this suppression is dependent on F. In 
particular, we see from Fig. 7 that the effect of disorder van¬ 
ishes when F —» 0, and a topological gap (albeit small) exists 
in that limit. This is consistent with our results in Secs. Ill A 
and IIIB on the scaling of A„ and with respect to F 
[Eq. (20) and Eq. (23), respectively]: since at small F, 
scales as F while r”' scales as F^, the pairing term eventually 
dominates over disorder at small F, producing a superconduct¬ 
ing gap. 

Figure 8 shows the dependence of the optimal values of 
tunnel coupling for which the topological gap on the wire is 



Figure 8. (Color online) The optimal value of F at which the topo¬ 
logical gap on the wire is maximum (Fop), as a function of disorder 


in the bulk SC 



The parameters of the wire are equal 


to those of Fig. 7(b), i.e., p = A, fi = 2A, and a = 



largest, as a function of the disorder in the bulk SC. We see 
that with a more disordered SC, it is actually more favorable 
to have a smaller coupling between the SM and the SC in 
order to generate a larger topological gap. We believe this 
nonintuitive finding (Fig. 8) to be an important result for the 
fabrication of optimal SC-SM hybrid structures for the real¬ 
ization of Majorana fermions— the tunnel coupling could be 
strong for an ultraclean SC (where / » ^), but for dirty SCs, 
one is far better off (as shown in Fig. 8) having a rather small 
SC-SM tunnel coupling. 


D. Discussion 


We have analyzed, within the framework of self-consistent 
Born approximation, the effect of disorder residing solely in 
the bulk of the SC on the spectral properties of the proximity- 
induced topological superconductivity in the SC-SM hybrid 
system. The dependence of the induced pairing gap and the 
induced disorder on the coupling strength is theoretically ex¬ 
plored. Crucially, we find that the topological gap induced 
on a SM wire with SOC and Zeeman splitting can be very 
susceptible to the disorder in the bulk SC when the SC-SM 
tunnel coupling is strong. While the specific optimal coupling 
strength depends on the details of the system, in general with 
high disorder in the SC a weak SM-SC coupling is preferable. 

These results have implications for the ongoing experimen¬ 
tal efforts to generate Majorana fermions by proximitizing a 
spin-orbit-coupled SM nanowire under magnetic field in con¬ 
tact with a SC. Although it is important to improve the inter¬ 
face quality between the two materials so as to generate a hard 
gap on the SM [58, 59], one should also be aware that a strong 
SC-SM tunnel coupling induces stronger disorder on the SM 
wire, if the SC is diffusive. Thus it is necessary to either use 
an ultraclean SC or to introduce a barrier between the SM and 
the SC so as to effectively reduce the coupling strength. 
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IV. CONCLUSION 


In this paper, we examined the effect of disorder in the bulk 
SC on a Shiba state or a proximate SM in the context of the 
current search for Majorana fermions in hybrid superconduct¬ 
ing systems. In both cases, we found that this type of dis¬ 
order can have signihcant detrimental impact, and could be 
an obstacle to create topological superconductivity in the hy¬ 
brid systems. In particular, disorder in the bulk SC can ran¬ 
domly shift the energy of the Shiba states in the ferromagnet- 
superconductor hybrid system, which is unfavorable for the 
existence of Majorana modes since realizing a Majorana- 
carrying topological system requires the fine tuning of the 
Shiba-state energy. (We mention that a complementary model 
[34, 70] of the ferromagnetic adatom chain on the supercon¬ 
ductor system, which is adiabatically connected [35, 73] to the 
Shiba model [69], assumes the system to be equivalent to the 
semiconductor nanowire on the superconductor structure, ex¬ 
cept for the spin splitting in the adatom chain being extremely 
large so that the system is completely spin polarized— in such 
a spin-polarized nanowire model of the ferromagnetic chain, 
the effect of bulk disorder is qualitatively similar in the semi¬ 
conductor and the ferromagnetic chain system with disorder 
being detrimental in the strong-coupling situation for reasons 
discussed in Sec. Ill above.) In the case of semiconductor- 
superconductor structure, the SM inherits both superconduct¬ 
ing pairing and disorder from the SC through the proximity 
effect. We find that the scaling exponent of inherited disorder 
with respect to the coupling strength between the two mate¬ 
rials is always larger than that of the inherited pairing. This 
implies that while the pairing term can dominate over disorder 
and produce a spectral gap on the SM at small coupling, upon 
increasing the coupling strength the inherited disorder will 


eventually dominate and destroy the induced SC gap. While 
the precise optimal value of the relevant coupling strength 
for producing the strongest topological superconductivity de¬ 
pends on the particular details of various parameter values, the 
key message of our theory for the choice of the most suitable 
topological materials parameters is that one should use ultra¬ 
clean bulk superconductors with extremely large normal-state 
low-temperature mean free path and tune the tunnel coupling 
to a suitable value lower than the bulk superconducting gap 
energy. We note that the disorder in the bulk superconduc¬ 
tor enters our theory through the dimensionless combination 

d - which implies that increasing either kpl or 1/^ 

in the parent superconductor should help to keep the disor¬ 
der effects weak in the system. This leads to our conclusion 
that among the commonly used parent superconductors in the 
experimental hybrid systems, probably Nb (Al) is the best 
(worst) choice, with Pb being somewhere in between since 
typically Nb (Al) has the shortest (longest) coherence length, 
thus making it easy (difficult) to satisfy the / » ^ condition. 
The detailed choice for the superconductor requires careful 
materials preparation with the longest (shortest) possible val¬ 
ues of the mean free path (coherence length) in the system. 
The precise prediction of our theory is simple: Choose a su¬ 
perconductor with the smallest possible value of the dimen¬ 
sionless disorder parameter “d” and, given this value of “d”, 
tune the tunnel coupling so that it equals Top, shown in our 
Fig. 8. 
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